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Abstract 

To determine the dissipation induced by the surface of fused silica fibers, 
we measured the quality factor of fibers having various diameters. We mea- 
sured a maximum quality factor of 21 million and extrapolated to obtain an 
intrinsic quality factor for fused silica of 30 million. Dissipation in the sur- 
face dominated at diameters less than about 1 mm. We developed a method 
for characterizing surface-induced dissipation that is independent of sample 
geometry or mode shape. 
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I. INTRODUCTION 



In experiments of sufficiently high sensitivity, fluctuations in the thermal energy of signal- 
band degrees of freedom can be an important source of noise. These fluctuations are due to 
coupling between the signal-band degrees of freedom of the detector and the thermal bath of 
other degrees of freedom. Since this coupling is also the cause of dissipation, one can relate 
the thermal noise in a detector to the dissipation of signal-band excitations. The relationship 
is quantified, for a wide class of couplings, by the Fluctuation-Dissipation Theorem: 



7T 2 f 2 
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where x(f) is the spectral density of the fluctuations, T is the temperature of the detector, 
and Z~ 1 {f) is the admittance of detector excitations. Si The real part of the admittance is 
proportional to the dissipation, so to minimize the thermal noise of a detector at a given 
temperature we must minimize its dissipation. 

In mechanical experiments, such as gravitational wave detectors, the relevant thermally 
excited degrees of freedom are actual mechanical vibrations of detector components. The 
level at which such a detector is limited by thermal noise can therefore depend on the 
internal friction of the material from which the detector components are fabricated. The 
lower the dissipation in the material, the lower the thermal noise. Fused silica (synthetic 
amorphous SiC^) has extremely low dissipation for audio frequency oscillations at room 
temperature. Therefore, fused silica is an excellent choice of material for the fabrication of 
critical components of a thermal noise limited detector, operating at room temperature in 
the audio frequency regime.! 

Many authors have reported measurements of the dissipation in fused silica.H However, 
the measured quality factors, Q, are usually limited by dissipation mechanisms other than 
the intrinsic dissipation of the bulk material. Particularly worrisome is the effect of 
a structurally defective and chemically impure surface, such as might result from normal 
handling, exposure to the atmosphere, or differential cooling during fabrication.!!! Although 
the current investigations are focused on fused silica, the condition of the surface is also a 
concern for other high-Q materials.0 

To investigate the amount of dissipation induced by the surface, we measured the dissipa- 
tion in fused silica fibers of varying diameters, and hence varying surface-to-volume ratios. 
In addition to providing information about the effect of the surface, our results may be 
extrapolated to give an approximate value for the dissipation in the bulk material. 



II. THE EXPERIMENT 

Our fibers were hand drawn in air, from Heraeus Co. "Suprasil 2" -brand fused silica rods, 
using a natural gas flame. The rods had diameters between 5 and 9 mm, and the fibers were 
drawn from the rods as supplied, with no additional surface preparation, either before or 
after drawing. 

The dissipation in the fibers was measured at room temperature using the ringdown 
method.0'0 Figure [IJ shows the experimental setup. The fibers, hanging freely, were excited 
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with a comb capacitor providing an oscillating electric field with a large gradient. Due 
to the dielectric properties of fused silica, an oscillating force is felt by the fiber. After a 
resonant mode was excited, the field was turned off, the capacitor grounded, and the fiber 
allowed to ring freely. 

The displacement of the fiber as a function of time was measured by a split photodiode 
shadow sensor. The envelope of the displacement was extracted and fit to the functional 
form 

Xoe -*fn<l>(fn)t + ^ ( 2 ) 

where t is time, xq is the initial amplitude, C is the level of sensor and amplifier noise, and 
4>{f n ) is the loss angle at the resonance frequency f n . (The quality factor of a resonance is 
equal to the reciprocal of the loss angle; Q = (plfn)^ 1 .) In this experiment the dominant 
source of noise was sensor and amplifier noise. Equation ^| does not include a term describing 
seismic excitation of the fiber, which unlike sensor and amplifier noise, adds to the signal in 
quadrature^ At the frequencies of interest , seismic excitation of fiber resonant modes was 
not detectable above the sensor and amplifier noise, and including such a term in Eq. did 
not measurably improve the quality of the fits. 

Great care was taken to eliminate all extrinsic sources of dissipation, or "excess loss." In 
any measurement of fiber Q's, there are several sources of excess loss that must be considered: 
residual gas damping, rubbing at the clamp-fiber interface, recoil damping, and eddy-current 
damping.!!^ 

The damping due to residual gas molecules in the vacuum chamber is given by 

4>gas = —ji (3) 

7TJ« P fiber 

where d is the fiber diameter, / is the frequency of vibration, v is the average speed of the 
gas molecules, p gas is the mass density of the gas, and p fiber is the mass density of the fiber. 
Gas damping was made negligible by conducting our measurements at pressures around 
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10 torr. Taking typical values for the parameters, we estimate 4> gas ~ 10 

To eliminate rubbing at the clamp-fiber interface we ensured that the fiber oscillation 
did not induce elastic deformation of the fiber material in the clamp. This was achieved by 
drawing the fiber monolithically from a much thicker rod (i.e. leaving the fiber attached to 
the rod from which it was drawn), then clamping the rod in a collet as shown in Fig. |l|.@ 

Recoil damping is due to coupling between the resonant modes of the fiber and low-Q 
resonances of the support structure. Resonant modes of the fiber having frequencies close to 
resonances of the supporting structure will be very strongly damped, while other modes will 
not be as strongly coupled and will have less recoil damping. The main diagnostic for recoil 
damping is therefore a strong frequency dependence of the Q.@ In this experiment, recoil 
damping was minimized by isolating the fiber resonances from the support resonances by a 
structure analogous to a double pendulum. This also serves to isolate the fiber from seismic 
excitation at the frequencies of interest. Figure shows the three different double pendulum 
type structures used in this experiment. The structures were monolithic to prevent interfacial 
rubbing, and in each case we measured the dissipation in the lowest fiber only. In the design 
of type 3, the lowest "fiber" is simply the undrawn rod, as supplied. 
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Eddy-current damping occurs when the oscillating fiber carries a charge, and the motion 
of the charges induces eddy currents in nearby conductors. Resistance in the conductors 
dissipates the mechanical energy stored in the currents, degrading the Q. The dissipation in 
any given fiber was not noticeably dependent on the arrangement of, or distance to, nearby 
conductors. We conclude that eddy-current damping was negligible in our measurements. 

Finally, we note that our fibers were drawn from fused silica rather than fused quartz 
rods. Fused quartz is fabricated from naturally occurring SiC>2, while fused silica is made of 
synthetic SiCV Fused quartz has been measured to have intrinsic room temperature Q's of 
at most a few million.0O'0'@ These values are too low for our purposes and would obscure 
the dissipation induced by the surface layer. 

Figure § shows the variations in Q measured with and without some of the precautions 
mentioned above. If the fiber is held directly in a clamp instead of being left attached to the 
rod from which it was drawn, friction due to rubbing in the clamp dominates. Without an 
isolation bob, recoil damping dominates and there is strong frequency dependence. Fibers 
clamped and hung by the rods from which they were drawn, and having a central isolation 
bob show the highest Q's, with much less frequency dependence. In these fibers, the differ- 
ence between fused quartz and fused silica becomes apparent, with fused silica exhibiting 
substantially higher Q. 



III. DISSIPATION VERSUS FREQUENCY FOR A TYPICAL FIBER 

For each fiber, the dissipation was measured at a number of resonance frequencies. For 
fibers of the type 1 design, the resonance frequencies agree well with the resonance frequen- 
cies of a beam of circular cross section clamped at one end: 



where Y is Young's modulus, d is the fiber diameter, L is the fiber length, p is the mass per 
unit volume, and n = 1, 2 ... is the mode number .0 Figure |] shows the agreement between 
the measured and predicted resonance frequencies for a typical fiber of type 1. The mode 
frequencies may be used to calculate the diameter of the fiber, and we find good agreement 
with the average diameter measured using a micrometer. Similarly, for the fibers of type 3, 
the resonance frequencies agree with those of a free beam of circular cross section. As 
expected, the resonance frequencies of the fiber of type 2 are not well modeled by either a 
free or a clamped beam. This is due to the large diameter of the excited fiber relative to the 
first (lowest) isolation bob which takes some part in the motion. Because of this, the fiber 
of type 2 was made with an extra isolation bob. 

Figure [| shows the dissipation versus frequency for the fiber whose resonance frequencies 
are shown in Fig|j. The graph shows a loss peak around 100 Hz which is due to thermoelastic 
damping.il Thermoelastic damping is given by 
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where Y is Young's modulus, a is the thermal expansion coefficient, To is the fiber temper- 
ature (absolute scale), C is the heat capacity per unit volume, and / is the frequency of 
vibrations. The constant a is 

2tt Cd 2 
13.55 k ' 

where d is the fiber diameter and k is the thermal conductivity. We found good agreement 
between the average measured diameter and the diameter calculated from the position of 
the thermoelastic damping peak. 

As can be seen from the graph, the dissipation is not purely thermoelastic and can be 
modeled by the form 

<PU) = 4>therrn(f) + 4>c, (6) 

where 4>therm(f) is thermoelastic damping and (fi c is a frequency- independent damping term. 
The dissipation in most of the fibers is in good overall agreement with the form of Eq. |6|. 
(Several fibers have a small number of modes showing anomalously high dissipation. This 
indicates an undiagnosed source of dissipation affecting these modes.) 

The appearance of a frequency-independent loss angle C suggests the presence of a 
dissipation source whose microscopic components have a wide range of activation energies.ii 
Dissipation in the bulk material is likely to be of this type. However, a frequency-independent 
term might also arise from defects or impurities in the surface layer, or possibly from other 
sources. 

One way of obtaining further information as to the source of <ft c is to measure its depen- 
dence on the fiber diameter. In this way it is possible to distinguish between dissipation 
occurring in the bulk volume of the fiber and dissipation occurring in the surface layer. 



IV. MODEL OF DIAMETER DEPENDENCE 

The constant loss angle C is modeled as consisting of two parts, one due to dissipation 
in the bulk and one due to dissipation in the surface layer: 

C = {AE mk + AE surf )/E, (7) 

where AE^uik is the energy lost per cycle in the bulk material, AE sur f is the energy lost per 
cycle in the surface layer, and E is the total energy stored in the oscillating fiber. 

If we make the rather general assumption that AE sur f is proportional to the surface area 
S while AEb u ik is proportional to the volume V, we may write 

AE sur f ^ S ^ 

AEbuik V 

The coefficient of proportionality depends on the ratio of the loss angle of the surface layer 
to the loss angle of the bulk material. It also depends on the ratio of energy stored in the 
surface layer to energy stored in the bulk. Some complication arises because both the loss 
angle of the surface layer and the density of energy stored in the surface layer are functions 
of depth. While we can normally calculate the energy density from the strain profile of the 
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mode shape, we know little about dissipation in the surface layer. Since we are interested 
in characterizing the dissipation in the surface independently of the mode of oscillation or 
sample geometry, we write the coefficient of proportionality as a product of two factors 

AE sur j ^ d s 

AE bu ik V I S 

The geometrical factor /x depends only on the geometry of the sample and on the mode of 
oscillation, while the "dissipation depth" d s depends only on the strength of dissipation in 
the surface layer relative to the bulk. The appropriate expressions for fi and d s are calculated 
in the Appendix. When the Young's modulus of the surface layer is the same as that of the 
bulk, we have 

1 r h 

d s = / (j)(n)dn, (10) 



4>bulk JO 

where n is a coordinate measuring the distance inward from the surface, <pbuik = ^E^uikl 1 E^uik 
is the loss angle of the bulk material, <p(n) is the loss angle of the surface layer as a function 
of depth, and h is the thickness of the surface. For our fibers, having circular cross section 
and oscillating in transverse modes, we have \x = 2. 

For samples with simple geometries, /i is of order unity and the volume-to-surface ratio 
has the same order of magnitude as the minimum thickness of the sample. When the 
dissipation depth is small compared to the minimum thickness of the sample, the effect of 
the surface on the dissipation is also small. When the dissipation depth is greater than or 
on the order of the minimum thickness of the sample, dissipation in the surface is likely to 
dominate. 

Since (f>(n) is seldom known explicitly, a measurement of the dissipation depth provides 
a convenient way of comparing the surface condition of different samples made of the same 
material. Since E « E^uik we may rewrite Eq. |7| in terms of the dissipation depth, 

<Pc = <Pbulk(l + ^^Jg)- t 11 ) 

In our case V/S = d/4, and the theory predicts 

<Pc = <Pbuik(l + 8^). (12) 

Equation [1^ is the model to which we shall compare our data on dissipation versus fiber 
diameter. 



V. DISSIPATION VERSUS DIAMETER 

For all of the fibers, the constant loss angle C was measured with the surface condition 
"as drawn". The surface of fibers drawn in a flame is largely free from microcracksili! 
and we tried to avoid damaging this surface. Although care was taken during transport 
and during installation in the apparatus, some of the fibers did get lightly knocked against 
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aluminum or glass components. This represents the only physical contact with the fiber 
surface. 

Measurements of <p c were made in the following way. Since </>(/) — > 4> c f ar from the 
thermoelastic damping peak, a measurement of the total loss angle at a frequency where 
thermoelastic damping is known to be negligible, constitutes a direct measurement of <p c . In 
each case <p c was taken as the lowest measured value of the loss angle for a particular fiber. 
In most cases, measurements of <fi could be taken at a sufficiently large range of frequencies 
so that those modes exhibiting the lowest dissipation gave a good approximation to the <p c 
asymptote. The three thickest fibers however, posed some problems. We were only able to 
excite two or three modes in each, and no direct verification could be made of the existence of 
a 0c asymptote. In addition, the dissipation in these fibers is very small, and correspondingly 
more sensitive to excess loss. Table | lists the Q's measured for the three thickest fibers. 
While Fiber G exhibits constant Q's, the Q's of Fiber J and Fiber K are quite frequency 
dependent. This is particularly striking for the split-frequency, third resonance mode of 
Fiber K, where the Q changes by a factor of 3 within 4 Hz. The source of this excess loss is 
undiagnosed. 

Figure ^| shows C versus diameter for the 10 fibers measured. Systematic errors are 
likely to be far larger than the uncertainty shown (which represents the repeatability). 
The main source of systematic error is the upward bias of the measured dissipation due to 
undiagnosed sources of excess loss. We fit the data to Eq. [12] and have tried to minimize 
the error induced by undiagnosed excess loss by including in the fit only those fibers whose 
graphs of dissipation versus resonance frequency do not have points deviating significantly 
from the form predicted by Eq. ^|. (For example, the fiber whose dissipation versus frequency 
data is shown in Fig. | satisfies this criterion well.) The fit determines d s and 4>b u ik, which 
have the values 

d s = 180 ± 20 /mi, (13) 
<p bulk = 3.3 ± 0.3 x 10~ 8 . (14) 

The relationship between d s and other measures of surface condition (fiber strength, surface 
roughness, impurities, etc.) is an interesting and open question. Further research is required 
to understand the dependence on surface preparation methods, storage times, manufactur- 
ing and handling, etc. The above value for 4>buik is consistent with the lowest dissipation 
measured in fused silica.!'! A quality factor of approximately 3 x 10 7 has been seen in 
hemispherical resonators of surface-treated, Russian brand KS4V fused silica at 3.7 kHz.!] 
With knowledge of the dissipation depth, we can use Eq. |], to calculate the fiber diame- 
ter, d eq , at which surface-induced dissipation becomes equal in importance to bulk-induced 
dissipation: 

d eq = 8d s ~ 1300 jum. 

In order to obtain an estimate for the average loss angle in the surface layer, we model 
the surface layer as a homogeneous shell of thickness h having a depth-independent loss 
angle, <p(r) = <p sur f- Equation |H] gives 

ds = /^W (15) 

4>bulk 
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The literature suggests several mechanisms for chemical surface damage penetrating to a 
depth of order 1 /im.0 Taking h = 1 /im and using the values given by Eqs. 13 and 14 we 
obtain 



■'surf 



10" 



(16) 
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APPENDIX: FORM OF THE GEOMETRICAL FACTOR AND THE 

DISSIPATION DEPTH 

It is possible (in a continuum approximation) to define a point loss angle, 

0(f) = Ap E (f)/p E (r), (Al) 

where f represents the location within the sample, p E {r) is the energy density stored at f, 
and Ap E {f) is the negative change in p E {r) per cycle. 

The bulk may then be defined as the region of points within the sample where 0(f) 
is constant. For the purposes of our model, we assume that there exists a surface layer of 
maximum thickness h where 0(f) varies, while elsewhere in the sample 0(f) has the constant 
value <f) bu ik. 

Recalling 

p E {r) = l -Y{r)e\r) (A2) 

where Y(r) is the Young's modulus and e(f) the strain amplitude, we can write 

AE surf = l -( cf){r)Y(r)e 2 (r) d 3 r (A3) 
2 Jsc 

where SC is the region constituting the surface layer. If we make the assumption that 0(f) 
and Y(r) may be treated as functions of depth alone, we can write 
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AE surf = I [ h <J)(n)Y(n) jjf ^ e 2 (r)rf 2 r| dn, (A4) 



2 Jo 



where n measures the distance in from the surface and S(n) is a surface of integration at 
depth n, parallel to the actual surface of the sample. For the bulk we have 

AE bulk = l<p bulk Y bulk J v e\f)d 3 r (AS) 

where V is the region of the bulk and Ybuik is Young's modulus of the bulk. 
We now make the assumption that h is sufficiently small that 

e(feSC)^e(feS), (A6) 

where S = <S(0) is the actual surface of the sample. Using Eqs. |A4| -|A6| we can write 

AE^-^V/S' (A7) 

where V is the volume of the sample, S is its surface area, and /i and d s are assigned the 
values 

V J s e 2 (r)d 2 r , A , 

b J v e 2 (r)d d r 



d s = - \z j h (j)(n)Y(n)dn. (A9) 

(Pbulk^bulk J 

The geometrical factor \x is a dimensionless constant which depends only on the sample 
geometry and on the class of resonances excited. The dissipation depth d s has the dimensions 
of length and provides a direct measure of the total dissipation induced by the surface layer 
(normalized to the dissipation in the bulk). In uncoated samples, the dissipation depth 
provides a measure of the physical and chemical damage suffered by the surface of the 
sample. 

For a transversely oscillating fiber of circular cross section, the strain amplitude is given, 
in cylindrical coordinates, by 

e (f)= ^\ J rcos6, (A10) 
oz 2 



where u(z) is the displacement amplitude of the fiber from equilibrium. Using Eq. [A^ we 
immediately have 

[i = 2. (All) 
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FIGURES 



FIG. 1. Schematic diagram of the experimental setup. The signal from the split photodiode 
shadow sensor is fed through a differential amplifier, bandpass filter, a lock-in amplifier and com- 
puter data acquisition system. 



FIG. 2. Three different designs for fiber isolation structures. In each case, only the dissipation 
in the lowest section was measured. In total, ten fibers were measured. Seven were of type 1, 
having diameters less than or equal to about 1050 /im, one was of type 2, having diameter of about 
3500 fj,m, and two were of type 3, having diameters of about 4890 /urn and 5930 /im. 



FIG. 3. The effect of isolation techniques. The left-hand graph shows the results for three 
natural fused quartz fibers having diameters between 350 and 500 p, The crosses indicate a fiber 
of type 1 (see Fig. |2|) clamped in a collet by the rod from which it was drawn. The circles indicate 
a similar fiber clamped in a collet by the rod from which it was drawn, but lacking an isolation 
bob. The triangles indicate a fiber detached from the rod from which it was drawn, lacking an 
isolation bob, and clamped between two plates. The right-hand graph shows the results for a single 
fused silica fiber of type 1, having diameter approximately 400 fim. The crosses represent the fiber 
clamped in a collet by the rod from which it was drawn. The circles represent the same fiber, in 
the same collet, but clamped and hung from the isolation bob. 



FIG. 4. Fit of theoretical cantilever beam resonance frequencies to the resonance frequencies 
of a typical fiber of type 1. The data shown is from a fiber of length 13 ± 1 cm and diameter 
120 ± 20 fim. (The relatively large uncertainties are due to the taper where the fiber exits the 
central bob.) The crosses are the measured values and the solid line is a fit of the measured values 
to Eq. § 



FIG. 5. Measured <j) vs. resonance frequency for a typical fiber of type 1. The circled bars 
represent the measured dissipation and uncertainty. The solid line represents the theoretical ther- 
moelastic plus constant damping and is a fit of the measured values to Eq. ||. The dashed line is 
the thermoelastic term from the fit. 



FIG. 6. Measured (fi c vs. average measured fiber diameter. The circled points represent fibers 
whose dissipation vs. frequency graph follows Eq. |6] without anomalous points. The uncertainty in 
(f> shown is the approximate repeatability (5%). The uncertainty shown in the average measured 
diameter is an estimate for the uncertainty induced by the fiber taper, where it exits the rod from 
which it was drawn. The solid line shows a least squares fit of Eq. 12 to the circled points. The 
diameter of the thinnest fiber was not measured directly but obtained from the position of the 
thermoelastic damping peak. 
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TABLES 

TABLE I. Quality factors exhibited by the the three thickest fibers. Repeated measurements 
typically vary by about 5%; this is larger than the measurement uncertainty. The large uncertainty 
in the average measured diameter is due to the taper at the upper end of the fiber. 
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Table Q 



Sample Avg. diameter (pm) Mode number / (Hz) Q 

Fiber G a 3500 ± 250 2 732.0 ±0.1 2.1 x 10 7 

3 1582.5 ±0.1 2.1 x 10 7 
Fiber J b 4885 ± 120 2 C 2157.3 ±0.1 1.0 x 10 7 

2 C 2167.3 ±0.1 0.41 x 10 7 

Fiber K b 5934 ± 70 3 C 1725 ± 0.5 0.68 x 10 7 

3 C 1729 ± 0.5 2.0 x 10 7 

4 3364.0 ±0.1 0.35 x 10 7 



a Type 2. 
b Type 3. 

c Since fiber cross sections are not perfectly circular, mode frequencies are split. 



14 



Figure 1 , Gretarsson et al., Rev. Sci. Instr. . 
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Figure 2, Gretarsson et al.. Rev. Sci. Instr. . 
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Figure 3, Gretarsson et al., Rev. Sci. Instr. 
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Figure 4, Gretarsson et al. ( Rev. Sci. Instr. 



4000 




2 4 6 8 10 12 14 16 

Mode number, n 



18 



Figure 5, Gretarsson et al., Rev. Sci. Instr. 




Figure 6, Gretarsson et al., Rev. Sci. Instr. 
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